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We derive the formula of the complex Doppler shift in a two-dimensional, dispersive metamaterial 
and we show that a moving, monochromatic radiation source generates multiple frequency modes. 
The role of the group velocity is stressed and the Doppler shifted radiation field exhibits features 
of the Cherenkov effect. The presented theory is also applicable to the case of a moving, non¬ 
oscillating charge and explains many peculiar characteristics of the Cherenkov radiation in left- 
handed metamaterials such as the backward direction of power emission, the constant radiation 
angle and the lack of velocity threshold. 


I. INTRODUCTION 

The negative refraction index media, first considered by Vesel^o [I|, and recently realized in the 
form of metamaterials @,9 exhibit many fascinating properties 9 9 potential applications such as 
imaging beyond the diffraction limit 9- such materials, commonly called left-handed media (LHM) 
due to left-handed triada of vectors E, H and k of propagating waves, many usual optical phenomena 
become reversed, however some of them are unreversed, just to mention the so-called rotational Doppler 
effect 9 as well as the spin Hall effect in LHM 9 • The class of reversed phenomena includes refraction 
9 9 1 Goos - Hanchen shift 9 and the Doppler effect [I| . The reversed Doppler shift has been observed 
in a variety of systems such as magnetic thin film 0 , photonic crystals [ll| , transmission lines [l^ and 
acoustical metamaterials [l3l| . In dispersive media, the so-called com plex Doppler effect occurs, when a 
monochromatic source generates wave modes of multiple frequencies [l4| . 

Another example of reversed phenomena in LHM is the Cherenkov effect [l|' This radiation, which 
was discovered experimentally in 1934 by Cherenkov [l5j | and later described theoretically by Frank and 
Tamm [l^ occurs when a charged particle moves at a speed exceeding the phase velocity of the waves 
in a medium and has been extensively applied for particle detection. In usual, right-handed materials, 
the angle of radiation is smaller than 7r/2, so that the energy is emitted in the forward direction. One 
of the important feature of negative index materials is a negative value of the phase velocity so the 
Cherenkov effect is reversed and the waves are emitted in backward direction [I|. It has been observed 
experimentally by Xi et al. m. The Cherenkov radiation in LHM has been studied by Lu et al. [l^ for 
both lossless and lossy cases and they found that maintaining a forward k vector of radiation in negative 
index media exhibits a backward emission. An elegant review of Cherenkov radiation in photonic crystals 
descri bing a backward-pointing radiation cone and backward direction of emission was presented by Luo 
et al. |l9| . More recently in [20| it was pointed out that the reversed Cherenkov radiation has a distinct 
advantage which allows the photons and charged particle to naturally separate in opposite directions, 
minimizing their interference. Nowadays there is a considerable interest in developing a new class of 
velocity-sensitive particle detectors based on left-handed materials [I^ which can be used to enhance 
the radiation [2lj or control its emission angle [^ . 

Here we present a simple, unified analytical description of the Doppler and Cherenkov phenomena in 
an idealized, dispersive medium. The description of the complex Doppler effect presented in our previous 
paper [l3 | is extended to the two-dimensional case. It is shown that a moving, monochromatic source 
generates two distinct frequency modes and the generated field exhibits features of both the Doppler and 
the Cherenkov effects. Starting from the first principles, we are able to explain the angular distribution 
of the energy for the fast and slow frequency modes of the two-dimensional, complex Doppler effect in 
LHM. Moreover, the reversed Cherenkov radiation emerges as a particular case of the presented theory 
as the source with the frequency wq = 0 is shown to be formally equivalent to a moving charge. The 
intrinsic feature of the reversed Cherenkov effect in the considered LHM is the lack of a source velocity 
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threshold above which this phenomenon occurs. This fact is explained in a general wa y a nd is shown to 
have a different origin than in the case of periodic media such as nanowire structures [23|. 

The theoretical findings are confirmed by numerical simulations based on the Finite Difference Time 
Domain (FDTD) method which is a common tool in the study of electromagnetic waye propagation 
in various media sm. The performed numerical simulations are in agreement with our theoretical 
predictions based on the first principle approach. 

This paper is organized as follows. In section 2 we present the theory of the complex Doppler effect in 
two-dimensional LHM, section 3 is devoted to the description of the simulation setup. A spatial structure 
of fields emitted by a moving source characterized by nonzero and zero frequency, describing the complex 
Doppler and the reversed Cherenkov effect respectively, is discussed in section 4. Conclusions are drawn 
in section 5. Finally, in Appendix A we prove that the Cherenkov cone angle varies only within small 
interval and is almost insensible to the source velocity. 


II. THEORY 

Consider an obseryer C receiving wavefronts emitted by a point source moving at a velocity v, where 
the initial angle between source-observer line and v is 0, as shown on the Fig. [T] The radiation source 
has a nominal frequency wq = I?- in its own reference frame. Assuming that the first wavefront was 



FIG. I: Source of radiation moving with velocity v emits two subsequent wave fronts at points A and 

B, towards observer C. 


emitted at t = 0, the time instant in which it will be registered is 

a + b n{uj){a + b) 

/ \ 5 

Vp[UJ) c 

where Vp(uj) is the phase velocity in the medium, and n(uj) is the refraction index. These quantities 
are observed in the reference frame of the detector and are functions of the detected frequency uj. The 
second wavefront will arrive at the receiver at the time instant 
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where (5 = v/c. The period measured by the observer will be 
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where the absolute value was used to ensure that the period is a positive quantity. In a typical case, the 
distance e traveled during a single period is negligible. Therefore, the angle a w 0 and b m d. By using 
the relation 

T 

a = ecos9 = v — , cos 0, (4) 

one obtains 
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similar to the result obtained in [^ . where in our case the absolute value is necessary to ensure that 
the frequency is a positive quantity. The above relation is implicit and has multiple solutions in a 
dispersive medium, leading to the so-called the complex Doppler effect [2^, where a monochromatic 
source generates multiple detected frequencies. 

To solve the Eq. [Sj one has to specify the dispersion relation of the medium. A remarkably simple 

description discussed by Veselago [Ij in the context of the left-handed media is given by the Drude model, 

^2 ^2 

with the permittivity and permeability given by ejoj) = 1 — and ^(tu) = 1 —respectively, where 
Wpe and ujprn depend on the metamaterial structure (27| . The proposed relations are applicable to a wide 
range of metamaterials iil- For further simplification, one can assume ujpe = ojpm = ^p, so that the 
dispersion relation for the refraction index n is 


n(w) = 1 - 


( 6 ) 


With the above dispersion model, from Eq. ([S]) one obtains two explicit solutions of the following form 
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where 


X = (3 cos 9. (8) 

The important property of the Eqs. © is the range of parameters where the frequencies oJi are real, 
indicating non-evanescent wave modes detectable in the far field of the source. By taking, for simplicity, 
Wq = ^ujp so that n{u}o) = —1, one obtains 


COS0 < 


2 - v'2 -k 2/32 
4/3 
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so that for sufficient source velocity, no propagating modes will be observed for the angle 0 smaller than 
some limit value 0 > 0o > tt/2. By inspecting the implicit Eq. ([5]), one can see that for the negative 
value of the refraction index, the Doppler shift is reversed, so that the frequency of waves emitted in the 
forward direction {9 < 7r/2) will be downshifted. In the considered model, the frequency of the upshifted 
modes at 0 > tt/2 is inherently limited; as a; ^ ujp, the refraction index n —>■ 0 and the denominator of 
Eq. ([S]) is equal to 1. Therefore, the frequency uip cannot be exceeded when wq < Wp. This means that 
a source having frequency wq such that n(wo) < 0 will generate only negative phase velocity modes with 
n(a;) < 0. 


III. SIMULATION SETUP 

The performed simulations were based on a standard FDTD algorithm [2^ where the Maxwell’s 
equations are solved in the time domain. A two-dimensional system with transverse magnetic field is 
chosen to facilitate the simulation of Cherenkov radiation the electric field has two components in 
the plane of propagation E = [Ex,Ey,0], and the magnetic field has a single component H = [0, 0, H^]. 
To model the dispersive medium, the Auxiliary Differential Equations (ADE) method is used [^. The 
implementation is a straightforward extension of the method used in our preceding work [l^ and is based 
on the time domain calculation of the polarization and magnetization in the medium. In the considered 
Drude model, these quantities are given by equations of motion P = uj^E and M = oJpH. 

The source with frequency ujq moving at a velocity V along x axis is modeled as a time varying current 
density in the form 


J = Jq exp 


— (* — Vt)'^ 
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The above representation describes an oscillating dipole oriented along x axis [T^ , where the usual Dirac 
delta function is replaced with sufficiently narrow Gaussian width, with the width a significantly smaller 
that the source wavelength A(a;o). The only relativistic effect taken into account is the time dilation 
affecting the frequency ojq, which is needed to keep the consistency with the relativistic description of 
the Doppler shift presented above. In the performed simulations, the Lorentz contraction and other 
factors affecting the field amplitude can be disregarded, as the following discussion of the radiation 
intensity is qualitative in nature and the relativistic effects at the considered speeds are not significant. 

The simulation space is divided into 300x300 grid and is surrounded by absorbing boundaries to 
reduce reflections. The frequency ojq was chosen to minimize the effect of the numerical dispersion and 
anisotropy [l^. This is especially important in the case of a dispersive medium described above, where 
the high frequency modes are limited by the finite time step, and the low frequency, short wavelength 
modes are constrained by the spatial step. 


IV. SPATIAL STRUCTURE OF THE FIELD 

Let’s assume that the moving source emits a quasi-monochromatic wave packet with central frequency 
Wo- According to Eq. ([7]), for any emission angle 0, there are two wave modes Wi(0), i = 1,2. Therefore, 
the wavefronts will propagate at an angle-dependent phase velocity Vpi(O) = Vpi[uj{6)\. However, their 
envelopes will move at a group velocity Vgi{9) = c/ (n(w) -f w|^) which is assumed to be positive. When 
Vg cos 9 < V, the emitted field will be detectable only behind the source, forming the so-called group 
cone [ 2 ^. In dispersive media, its apex angle may differ significantly from the usually considered wave 
cone formed by the wavefronts. The situation is shown on the Fig. [7] Assuming that the medium is 
isotropic and lossless, the phase velocity will be parallel to the group velocity. As it follows from the 



FIG. 2; Radiation cone formed by a single mode emitted at a given angle 9. Dashed line indicates the 

wavefront. 


geometry of the Fig. [2l the cone half apex angle a can be expressed as 


tan a = 


Vg (9) sin 9 
V — Vg{9) cos 9' 


( 11 ) 


To deduce the whole radiation pattern of the moving source, one has to consider the wave modes for all 
angles 9. An elegant way of presentation was pointed out in [I^ and it consists of plotting the group 
velocity Vg{9) as a contour and comparing it to the source velocity. By doing this, one can visualize the 
distance traveled by the source and the wave envelope at some given time, recreating the spatial structure 
of the radiated field. An example of such construction for /3 = 0.3c and n(wo) = —1 is shown on the Fig. 
[3a] and the corresponding field snapshot is presented on the Fig. [SB The solutions of Eq. © define two 
contours Vgi(9). The high frequency modes wi(0) are comparable to the nominal frequency wq and are 
characterized by relatively high group velocity Vgi{9). This solution can be associated with the usual 
Doppler effect; the backwards moving waves {9 > 7r/2) are upshifted, with their group velocity exceeding 
the value of V^(a;o) = c/3. Also, the transverse Doppler shift for the waves emitted in a direction 
perpendicular to the source motion is caused only by the relativistic effects, so in the low velocity regime 
lui{tt/2) ss ujq. For the forward moving waves, the frequency and the group velocity drops quickly as 
9 ^ 9q, possibly becoming much smaller than V and leading to the Cherenkov-like radiation pattern. 
The second wave mode UI 2 is characterized by small group velocity. A characteristic point of this solution 
is 0 = 7r/2, where W 2 = 0. As the angle approaches this value, the frequency becomes arbitrarily small. 
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FIG. 3: (a) Radiation cones of the two Doppler modes at /? = 0.3c and n{uJo) = —1. The angles 9i and 
62 mark the direction of radiation forming the cone, (b) The group velocity contour imposed on the 

field snapshot calculated with FDTD method. 


In a realistic scenario, the amount of the downshift is limited to the range where the dispersion model 
in Eq. (jO]) is applicable. Moreover, the low frequency modes are more affected by the absorption [l3 |. 
making their detection difficult. 

The space taken by the radiated field is limited by the maximum value of cone angle a measured 
between the source velocity and a line tangent to the group velocity contour (Fig. [2a|). Its value can be 
determined by solving diana/dO = 0, which leads to 

dV 

If = sin e + VV„ cos 9. (12) 

y 39 

The angle 9 obtained from the above equation marks the point of intersection and, according to the 
Fig. m describes the direction of radiation forming the cone. As it was shown in [ 2 ^ at this angle, the 
maximum of the radiation maximum intensity is expected. This is due to the fact that for all angles 
near 0, the cone angle a is almost the same, so that constructive interference occurs. In a dispersionless 
medium (e. g. Vg = Vp and dVgjd9 = 0), the Eq. (1121) reduces to the classic relation for the Cherenkov 
radiation angle 

cos 9 = ■^. (13) 

In the case of a dispersive material, the Eq. m can be solved separately for the two modes Vgi{9), 
yielding characteristic angles 9i. These critical angles, as well as and the associated cones described by 
Eq. (HU are shown on the Fig. [3U One can see that the first, wide cone is formed by the high frequency 
modes wi. The radiation is emitted at an angle 9i < 7r/2 and the wavefronts are roughly parallel to the 
cone surface, in a manner similar to the Cherenkov radiation (Fig. I3bl) . The low frequency modes UJ 2 
form a much narrower cone. The emission angle 02 > 7r/2 and the wavefronts are almost perpendicular 
to the cone surface which is a characteristic feature of the reversed Cherenkov effect [^. The first 
contour Vgiifi) becomes almost circular for 9 > 9i. In this area, the radiated field is very similar to the 
case of the Doppler effect in dispersionless medium - a detector positioned at an angle 9 will register 
only single wave mode, so that the wave frequency lu(9) can be readily measured. On the other hand, as 
the angle approaches 9o, the frequency of the both modes changes quickly, so the field measured along 
some finite angle range around 9q has a wide spectrum. Further analysis of this case is presented on the 
Fig. m The field snapshot from the FDTD simulation is presented in Fig. Clearly, the obtained 
cones divide the radiated field into three distinct regions. The highest field amplitude is generated by 
the modes traveling along the direction 9i and forming the outer cone. The inner cone is formed by the 
1 ^ 2 ( 02 ) modes. Another significant direction is 9q, where the field is a superposition of many modes with 
a wide range of group velocities, so that it is spreading over time. These findings are confirmed by the 
calculated time-averaged value of the Poynting vector S = E x H shown on the Fig. I4bl One can see 
that the initial radiation pattern shows a wide peak at 0 ss 0i which is formed by the fast modes uii. 
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After the time needed for these modes to leave the simulation area, the intensity peaks formed by the 
slow modes uj 2 can be observed. They are centered around Oq and 62 - As it is shown on the Fig. 03 
there is significant difference in the absolute value of the phase velocity between the fast modes wi and 
slow modes UJ 2 - In the latter case, it is significantly smaller than V. Finally, the measured frequencies 
presented in Fig. I4dl are also in a good agreement with the predicted value. To illustrate the influence 
of the source velocity on the field structure, another simulation was performed for (3 = 0.1. The results 
are shown on Figs. 03 and 133 In this case, the group velocity of the first mode wi is always significantly 
higher than V, so that no cone is formed and the radiation pattern differs only slightly from the field of 
radiating, stationary dipole. The second mode UI 2 is characterized by a very low group velocity, forming 
a very narrow cone, barely visible in the field snapshot in Fig. 03 

It is important to note the difference between the obtained radiation patterns and the Cherenkov effect. 
Here, the source of radiation is characterized by the nominal frequency wq. Moreover, the absolute value 
of the phase velocity of the first mode is significantly larger than V (Fig. 03, so that the condition in Eq. 
(03 cannot be met. One can also see that the cone angle 9i < 7r/2 which characterizes a non-reversed 
Cherenkov effect, despite the fact that the phase velocity is negative. Finally, for 0 > 0i, the cone 
transforms in a continuous manner into the typical Doppler shifted field pattern, generated by radiating 
point source. On the other hand, the second, narrower cone is formed by the modes characterized by 
1 1 ^ 2 ! < V and Vg 2 < V; the radiation pattern is simpler, closely resembling the case of the reversed 
Cherenkov effect. 

To compare the results to the generic Cherenkov radiation, our source model can be adapted to 
simulate a moving charge by assuming cuq = 0. Under such condition, the Eq. m describes wave modes 
with frequency uj = k ■ V = nfi cos 0, generating the spectrum of a point charge moving with velocity V. 
A particle of small but finite dimension moving along x axis can be modeled as a charge density [32j| 


p{x, y, t) = po exp 


— (x — Vt)'^ 

2 ^^ 


exp 



which corresponds to the source current [if 


J = PqV exp 


— (x — Vt)'^ 
2^2 


exp 



(14) 


(15) 


As a test case, a simulation was performed in a dispersionless dielectric characterized by n = 4. The 
charge moving at a velocity U = 0.3 c is expected to create a cone at the angle 0i = arccos(I^/U) ~ 34°. 
The geometrical construction based on the group velocity contour is presented on the Fig. [Sa] The 
field snapshot is shown on the Fig. I5bl The Cherenkov cone is clearly visible and its angle matches the 
theoretical prediction. In the case of a dispersive medium described by Eq. m, the Doppler shift Eqs. 
© remain valid for wq = 0, reducing to a single solution 


One can check that for any angle 0, the phase velocities of the wave modes described by above relation 
fulfill the Eq. (fUTll . This important property should be stressed; in the classical Cherenkov radiation, the 
single radiation angle is a function of the two constants - the phase velocity and the source velocity. Here, 
the phase velocity of the emitted mode varying, being a function of the source speed and angle 0. In 
other words, the radiation pattern is a superposition of cones corresponding to different group velocities 
0. The real value of the frequency associated with propagating waves is obtained from Eq. (1161) only 
for 0 > 00 = 7r/2. This means that the Cherenkov effect is reversed which is consistent with the negative 
value of the refraction index [ij. The group velocity contour for V = 0.3c is shown on the Eig. [53 In 
contrast to the Cherenkov effect in dispersionless medium, the frequency range of the emitted modes is 
limited, and their group velocities vary from U, = 0 at 0 = 90° to I4 ~ 0.19 c at 0 = 180°. Due to the 
low values of the group velocity, the outer cone is relatively narrow, and the radiation angle 0i « 120°. 
Again, the structure of the field in the snapshot presented in Fig. I5dl is in an excellent agreement with 
the analytical prediction. As expected, the radiation pattern presented in Fig. 03 shows a narrow initial 
peak at 0 = 01, followed by the slow modes at 0o > 0 > 0i. When a single frequency is measured, the 
radiation angle agrees with the generic Cherenkov condition given by Eq. (1131) for the appropriate value 
of n(w). Therefore, the performed simulations show a good agreement with prior simulations where a 
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(a) 


(b) 





FIG. 4: Simulation results for the complex Doppler effect at ^ = 0.3 and n(a;o) = —1. (a) Field 
snapshot at t=500. The cones clearly divide the space into distinct areas, (b) Angular distribution of 
energy for the fast modes centered around and the low frequency modes at 0o and Qi. (c) The phase 
velocity contour. The radius of the black circle is equal to the source speed, (d) The frequency as a 
function of angle measured in the area where the field is monochromatic, (e) Field snapshot for 
/3 = 0.1. (f) The group velocity contour for (3 = 0.1. 


window function has been used to reduce frequency range [2^ and with experimental data obtained by 
Xi et aZ.[l^. 

Surprisingly, the dependence of the cone angle on the source velocity is insignificant; the result for 
V = 0.8c shown on the Fig. [Sjis almost indistinguishable from the case of D = 0.3c. This is confirmed on 
the Fig. [5^ where the numerically computed relation 0i(/3) based on Eq. (fT^ is shown. This interesting 
property depends on the dispersion relation of the medium and has been observed in various physical 
systems; for example, in the so-called Kelvin wake in water waves [^ . the group velocity is directly 
proportional to the phase velocity, which is given by Eq. m- This causes the Eq. (Ha to independent 
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FIG. 5: Cherenkov radiation for f3 = 0.3. (a) Group velocity contour in dispersionless dielectric with 
n = 4. (b) Field snapshot at t=500, showing typical Cherenkov cone, (c) Group velocity contour in 
left-handed medium, (d) Field snapshot at t=500, showing reversed Cherenkov effect, (e) Radiation 
pattern for /3 = 0.3, exhibiting a narrow peak at Oi. (f) Radiation pattern for /3 = 0.8. 


of the source speed; the details are presented in Appendix A. 

One can see another significant departure from the typical Cherenkov effect in a dispersionless medium; 
the Eq. (fTHll has a solution for any 1 > /? > 0. Therefore, there is no source velocity threshold above which 
the phenomenon occurs. Unless additional constraints are put on the dispersion model, the frequency 
and phase velocity of the wave modes can become arbitrarily small, matching the source velocity in 
accordance with the Eq. (11311 and producing radiation at any /3. Such radiation without a velocit; 
threshold has been observed with metallic gratings [3l| and predicted to occur in periodic structures 
such as nanowire metamaterial [2^ . The Cherenkov radiation frequency varies considerably with source 
velocity (Eig. I6bl) but, in general, it is comparable to the resonant frequency of the medium Wp, so that 
it is defined by the material properties 321. 























FIG. 6: (a) The dependence of the cone angle on the source velocity, (b) The dependence of the wave 

frequency on the source velocity for 9 = 120°. 


V. CONCLUSION 

We presented an unified analytical description of the two-dimensonal Doppler and Cherenkov effects 
in an idealized, dispersive metamaterial. It was shown that a moving, monochromatic source generates 
two distinct frequency modes and the generated field exhibits features of both, the Doppler and the 
Cherenkov effects. We prove, starting from the first principles, that the reversed Cherenkov radiation 
emerges as a particular case of the presented theory. The characteristic peculiarities of the reversed 
Cherenkov effect in LHM, namely the backward direction of emission, the constant angle of maximum 
radiation intensity and the lack of a source velocity threshold above which this phenomenon occurs 
are explained on the basis of presented theory and confirmed by numerical simulations. Our results 
are in agreement with measurements obtained in the experiment performed by Xi et al. 0. All these 
interesting features of the Cherenkov radiation in metamaterials characterized by Drude-like dispersion 
model might possibly provide a new way of frequency-based charged particle speed measurement, with 
the mechanical simplicity enabled by the constant radiation angle. 
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Appendix A: The angle of Cherenkov radiation 


In the Drude model medium described by Eq. m, the refraction indices for the phase and group 
velocities are given by 


n = l -1 ng = l + -|, 






(Al) 


therefore, these velocities are linked by the relation 




V COS0 


^9 2 ^ — 1 2/3 cos 0 — 1’ 


(A2) 


where the phase velocity is defined by the generic Cherenkov radiation condition in Eq. dT51) . The Eq. (fT^ 
defining the angle of maximum radiation intensity can be simplified to the following form 


cos^ 0(3 — 2/3 cos 0) — 1 = 0. 


(A3) 


Considering the fact that the real frequency modes are given by Eq. m only at angles bigger than 90°, 
for any given source speed P, the above relation has one proper solution at 0 « 120° (Fig. [7]). One can 
check that in the limit of /3 ^ 0, 0 w 125° and for ,0 = 1,0 = 120°. 



FIG. 7: Plot of Eq. I A3 1 for selected values of /3. 


It is interesting to mention that in the case of the water waves, one has 

1 V^cos0 
^^ 9 - 2^P- 2 ■ 

With such defined group velocity, the Eq. m takes a simple form 

sin^ 0 = 


with a solution describing the well-known Kelvin wake angle 0 ~ 35.26°. 


(A4) 


(A5) 
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